Asymmetry in self-assembled quantum dot-molecules made of identical InAs/GaAs 

quantum dots 
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We show that a diatomic dot molecule made of two identical, vertically stacked, strained 
InAs/GaAs self-assembled dots exhibits an asymmetry in its single-particle and may-particle wave- 
functions. The single-particle wave function is asymmetric due to the inhomogeneous strain, while 
the asymmetry of the many-particle wavefunctions is caused by the correlation induced localization: 
the lowest singlet ^Eg and triplet '^E states show that the two electrons are each localized on different 
dots within the molecule, for the next singlet states ^^E^ both electrons are localized on the same 
(bottom) dot for interdot separation d > 8 nm. The singlet-triplet splitting is found to be ~ 0.1 
meV at inter-dot separation d=9 nm and as large as 100 meV for d=4 nm, orders of magnitude 
larger than the few meV found in the large (50 - 100 nm) electrostatically confined dots. 
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PACS numbers: 73.63.Kv, 85.35.-p, 73.23.Hk 

Quantum dot-molecules (QDM) occupied by two elec- 
tronic spins have been proposed as a basis for quantum 
computation.^ Loss and DiVincenzc^ proposed a "swap 
gate" based on a simplified model in which the two lo- 
calized spins have Heisenberg coupling H — Js-tSi ■ S2, 
where Si and 5*2 are the spin-i operators for the two 
localized electrons, and Js-t is the effective Heisen- 
berg exchange splitting, being the difference in energy 
between the spin-triplet state with total spin S = 1 
and the spin-singlet with S* = 0. Successful operation 
would require a large-singlet triplet splitting Js-t (fast 
swap time^), and that the probability Qtot of the two 
electrons in state v simultaneously occupying one dot 
be small (maximizing entanglemenli^) . The search for 

a nanosystem with large Js-t and small Q^^t involves 
engineering of the properties of the corresponding many- 
particle wave functions. In a simplified molecular model, 
the single-particle wave functions of the two electrons 
are given by bonding Tpa^ = (xt + Xb)/^/^ and anti- 
bonding V'(T„ — {xt — Xs)/v^ molecular orbitals, con- 
structed from the individual orbitals x of the top (T) 
and bottom (B) dot. The many-particle states are then 
the corresponding product states Icr^o"^) ^ (singlet), 
and \cr^gaX) ~ (triplet), etc. In this picture both 
single-particle molecular orbitals and the many-particle 
states are delocalized on both dots constituting the dot- 
molecule. Here, we discuss via atomistic single-particle 
and many-body calculations two important deviations 
from this simplified molecular picture, leading to asym- 
metries both in the single-particle molecular orbitals due 
to the inhomogeneous strains, and in the many-body 
states (i.e. localization either on T or on B) as a conse- 
quence of correlation. For the many-body states we find 
Mott-like transitions for the first and third singlet states: 
both electrons are localized on one dot at large d and 
delocalized over both dots are small d. The double occu- 
pancy qIoj of the first singlet state is surprisingly large 
(~40%) for an interdot separation of 5 nm. The triplet 
states and the second singlet state are Mott-localized at 



every interdot separations with Qtot— (^o double occu- 
pation) for the triplet states and Q^^J^=1 for the singlet 
state. 

Previous models of dot-molecules have focused on 
electrostatically-confined dots^*^, which have a very 
large confining dimension of 50 - 100 nm. Such dots 
exhibit typical single-particle level spacings of Ate =3 - 
5 meV, Coulomb energies Jge of about 5 meV > Asg, 
whereas exchange energies and correlation energies are 
around 1 meV. Many experiments were recently done 
on two such coupled dots^'®, showing that the splitting 
between the bonding and anti-bonding molecular lev- 
els is as large as 3.5 meV at an interdot separation of 
2.5 nm, comparable to the single-particle energy spac- 
ing Aeg of the single dot. Because of the very large 
size of such dots, their single-particle levels can be de- 
scribed by simple one-band effective mass "particle in 
a box" models using the external potential generated 
by a combination of band offset, the gate potential and 
the ionized impurities&iS Alternatively, one can sim- 
ply assume a particle-in-a-parabolic well modelpii*i^ In 
these descriptions, the single-particle states are sym- 
metric , but many-electron symmetry breaking is pos- 
sible due to correlation effects, as shown via unrestricted 
Hartree-Fock treatment of the effective-mass approxi- 
mation (UHF-EMA)fi^, configuration- interaction treat- 
ment of the effective-mass approximation (CI-EMA),— 
or Mott-Hubbard model>i^ 

Here, we discuss localization effects and singlet- 
triplet splitting of electrons in vertically coupled self- 
assembled InAs/GaAs quantum dot molecules grown 
epitaxially^^'-'^^-ii. Such dots have much smaller confin- 
ing dimensions (height of 3 - 5 nm), and when made of 
InAs/GaAs their electronic level splitting is Aeg ~ 40 - 
50 meV, larger than the interelectronic Coulomb repul- 
sion Jee ~ 10 -20 meV, and the exchange energy K^^, ^ 2 
- 5 meV. In this paper, we show that in vertically aligned 
self-assembled dots, one can achieve singlet-triplet split- 
tings of up to 100 meV. 

Figure QJl shows the geometry selected for the dot 
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(a) Hydrostatic Strain 




Wetting layer 



(b) Single-particle electron level of QDM 
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(c) Molecular-orbital Model 
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FIG. 1: (a) (Color online) Contour plot of the hydrostatic 
strain Tr(e) in the two vertically coupled quantum dots. The 
inter-dot distance d is measured from one wetting layer to the 
next. The iso-strain values are also marked in the figure, (b) 
Molecular orbital energy levels vs. inter-dot distance. The 
dashed line is the average of ag and ct„. Ssp is the single dot 
s-p energy level splitting, (c) Sketch of bonding-antibonding 
splitting. St, Sb are "s" while Pt and Pb are "p" single dot 
orbitals on top and bottom dots. Here St + Sb ~ o-g and 
Pt + Pb ~ ttu are bonding states, while the St — Sb = o"u 
and Pt — Pb = Tfg are antibonding states. 



molecules, consisting of two-dimensional InAs wetting 
layers, a pair of 2 nm tall InAs dots in the shape of trun- 
cated cones embedded in a GaAs matrix. The considered 
dots are identical and describe the experimental extreme 
case of perfect growth. In reality, both dots can be ge- 
ometrically and compositionally different. This extrin- 
sic inequivalence would reinforce the intrinsic asymme- 
try of the system which our results describe. The hydro- 
static strain field Tr(e), with the iso-strain lines shown in 
Fig. ^i, is calculated atomistically by relaxing the bond- 
stretching and bond-bending forces according to the va- 
lence force field model (VFF).^*'^^ It is clearly seen in 
Fig. ^ that both dots have large and nearly constant 
hydrostatic strain inside the dots which decays rapidly 
outside the dots. However, even though the dots com- 
prising the molecule are the same, the strain on the two 
dots is different since the molecule lacks inversion sym- 
metry. We see that the top dot is slightly more strained 
than the bottom dot. Furthermore, the GaAs region be- 
tween the two dots is much more strained than in other 
parts of the matrix. 

Having established a realistic geometry and the re- 
laxed atomic positions {Rm^a}, we calculate the single- 
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(b) CI calculation of 
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FIG. 2: (a) The effective single-particle energy levels ex and 
bb of dot-centered orbitals on the top dot and bottom dot 
respectively, (b) Energy of two-electron states calculated from 
CI using all confined molecular orbitals. 



particle electronic structure by constructing a pseudopo- 
tential V{r) — J2m. a '^ai'"^ ~ '^■m,a) from a superposition 
of screened atomic potentials Va of species a= Ga, In, 
As. Here, Va is constructedSS by fitting to available ex- 
perimental data the bulk InAs and GaAs band energies, 
effective-masses, hydrostatic and biaxial deformation po- 
tentials, and band-offsets. The pseudopotentials used in 
present work are taken from Ref. 12^. The Schrodinger 
equation is solved by the Linear Combination of Bulk 
Bands (LCBB)2i method in a basis {^^f^ «(r)} of Bloch 

orbitals, i.e. ^i{v) = En,k,A C'i!k,A '/'i^k,'?('^) of band in- 
dex n and wave vector k of material A (= InAs, GaAs). 
The basis functions are strained uniformly to constant 
but different strains e*. We use *e = for the (unstrained) 
GaAs matrix material, and an average *£* value from VFF 
for the strained dot material (InAs). For the InAs/GaAs 
system, we use n = 2 for electron states on a 6x6x28 
k-mesh. Unlike the effective-mass description of single- 
particle state used for electrostatic dot^iiSiiiii^ here we 
allow inter-band and inter-valley coupling. 

Asymmetry of single-particle states: Figure ^ shows 
the single-particle dot-molecule energy levels as a func- 
tion of inter-dot distance d. These results can be gen- 
erally understood using bonding/antibonding molecular 
orbitals as shown in FigCJ;. At large d, the energy lev- 
els converge to those of single dots, showing an s-level 
and, at Ssp — 106 meV higher, two p-levels split by a 
few meV reflecting the atomistic C2v symmetry of the 
cylindrically symmetric dots^Si. Bonding and antibond- 
ing molecular orbitals form when the two dots interact: 
the two single-dot s orbitals form molecular ag and ct„ 
orbitals, whereas the four single-dot p orbitals (two on 
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each dot) form two molecular 7r„ and two molecular iTg 
orbitalspi^i Note that the splitting of ag from cr„ is not 
symmetric, as can be seen by looking at the average of 
both energies in Fig. ^ (dashed line). The average en- 
ergy increases with decreasing interdot separation in re- 
sponse to the strain exerted by the presence of the other 
dot (Fig^). Beyond this overall effect, the strain field is 
different on both geometrically identical dots because the 
dots are non-spherical (Fig^) and lack inversion symme- 
try. This causes an asymmetry of the molecular orbitals 
at short interdot distances. At d ^ 3.4 nm, we found that 
the bonding (anti-bonding) states are slightly more local- 
ized on the bottom (top) dot. In previous effective-mass 
calculations strain effects were not included 
and single-particle asymmetry was not found. Partic- 
ularly, in Refj^, an ad hoc parameter was introduced to 
force the asymmetry of the dot molecule. 

A more quantitative analysis of the asymmetry can be 
given when the molecular orbitals are transformed via a 
Wannier-like transformation into single dot states. These 
can be obtained from a unitary rotation of molecular or- 
bitals ipi, 

i 

where, ipi is the z-the molecular orbital and xi,p £^re the 
rotated, dot-centered orbitals (the l-th orbital localized 
on p= T or B dot). U are unitary matrices, U^U = I, 
chosen to maximize the total self- Coulomb energy 
Once U^'l are known, we define the orbital energies of 
the dot-centered states xi,p ^s: 

e/.p = {Xi,p\f\xi,p) = Yj^uly ' (2) 

i 

where T is the kinetic energy operator, and is the en- 
ergy of the i-th molecular orbital. The energies ei^p of the 
dot-centered orbitals are depicted in Fig.|3i. as a function 
of the interdot separation. We see that the single-particle 
energies for both B and T orbitals rise quickly as the 
inter-dot distance is reduced, but the energy of the top 
dot orbital raises faster. At d ~ 3.4 nm, there is an en- 
ergy splitting of ^ 12 meV between top and bottom dot 
orbitals, which causes the asymmetry of the wave func- 
tions between top and bottom dots. 

Many-particle symmetry-breaking: Having obtained 
the "molecular", single-particle energy (Fig. ^), and 
wave functions, we calculate all interelectronic Coulomb 
and exchange integrals J and K of ipi by numeri- 
cal integration^^, and set up a screened configuration- 
interaction expansion^. A microscopic position- 
dependent dielectric screening^ is applied to both 
Coulomb and exchange integrals to represent the inner 
electrons that are not calculated explicity. Considering 
six molecular orbitals ag, cr„, 7r„, iVg of Fig.^, we have 
a total of 66 Slater determinants. The many-body wave 
functions ^'j/ are written as linear combinations of these 



determinants |$c) as, ^fi, = J^c'^i^i'^) l*c)- The result- 
ing many-particle energies are shown as a function of in- 
terdot separation in Fig.l^b- The energy splittings Js-t 
between the ground state singlet ^S^"-* and triplet 
ranges from - 100 meV and is much larger than in elec- 
trostatic dot molecules (< 1 meV)i^i^ In Fig|3Ji we de- 
compose the two-electron wave functions into the leading 
configurations $i = Wjal), $2 W^f^l), ^"3 = Wj^^): 
and $4 — Icier J;). The ground state is the singlet 

^Sg"' state, followed by the three- fold degenerated triplet 
states -^S (we depict only the Sz=0 state made of $1 -I- $2 
in Fig|2J) and the next singlets -^E„ (made of $1 — $2) 
and ^Eg''-'. To explore the localization of these states, 
we plot in Fig|3|3 and|3|;, the pair correlation functions 
Pi,(ro,r) = |^'i,((ro, r)p where Tq is fixed at the cen- 
ter of the bottom dot. Pi,(ro,r) gives the probability of 
finding the second electron at position r given that the 
first electron has been found at Tq. For the ground state 
singlet ^T,g°'\ we see that at the small interdot separa- 
tion d—A nm, the probability to find the second electron 
in the top or the bottom dot are comparable, suggest- 
ing a molecular-like delocalized state. Accordingly, the 
wave function analysis reveals a dominant contribution 
from the product of two delocalized molecular orbitals 
$3. With increasing interdot separation, the electrons 
show correlation induced (i.e. the coupling between $3 
and <I>4) localization. At d~7 nm, the second electron 
is almost entirely localized on the top dot as shown in 
Fig. Et. A similar delocalized to localized transition ap- 
plies for the ^Eg''-* state, with the difference, that at large 
d both electrons are localized on the same (bottom) dot. 
In contrast, the triplet states ''E and the singlet ^E„ 
show localization at all interdot distances. However, for 
the triplet states the two electrons are localized on dif- 
ferent dots while for ^E„ both electrons are localized on 
the same dot. 

To understand the (hidden) broken-symmetry (i.e. the 
wave function localization) in the many-particle CI wave 
functions, and to study the degree of localization quan- 
titatively, we resort to the Wannier-like transformation 
of Eq. . The CI matrix elements expressed initially in 
the molecular basis tfji [Eq. ^] are transformed into the 
dot-centered Wannier basis xi.p [Eq. Q]- For example, 
ilx'ipiX'i'p')} denotes the configuration where one elec- 
tron is on the ^-th orbital of the p dot with spin a, and 
the other electron is on the I'-th orbital of the p' dot with 
spin a' . The two electrons can be either both on the top 
dots, or both on the bottom dots, or one on the top and 
the other on the bottom dots. We can thus define a "bi- 
electron localization parameter" Q^p as the probability 
of two electrons occupying the dot p— (T or B) at the 
same time in the many-particle state i^, 

Qip^= E P'^i\Xi,p,xC))' (3) 

la, I' a' 

where Pu{C) is the weight of the configuration C in the 
many-body wave functions of state ly. The total prob- 
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FIG. 3; (Color online) Panel (a) gives the weights of con- 
figurations <l>i, "I>2, "I>3 and <1>4 of the many-particle CI wave 
functions. Panels (b) and (c) depict the probability of finding 
the second electron at position r given that the first electron 
has been found on the center of the bottom dot (indicated by 
the arrows) for the interdot distances (b) d=4 nm and (c) d= 
7 nm. 



ability of two electrons being on the same dot is then 
Qtot — Qtt + Qbb fo'^ ^"th state. Figure^ shows the 
bielectron localization parameter Qp^J of Eq.Q for the 
many-particle states v— ^Eg, We see that: (i) For 

the ground state ^Eg"^ (Fig-S^), Qbb > Qtt at all inter- 
dot separations (as is the case for the state), whereas 



for state ^Sg ' (Fig.^), we always have Qtt > Qbb, in- 
dicating symmetry breaking for these many-body wave 



(a) 



has a very small 
nm) , whereas for 



functions, (ii) The ground state ^E 
Qpp at large inter-dot separation {d > 

^Eg*"-*, the probability of two electrons being on the same 
dot is close to 1. At smaller d, Qtot increases rapidly for 



Eg"'', while it decreases for "'^Eg'''*. (iii) Our calculations 
show that the coupling between ■^E and higher triplet 
configurations is negligible, thus the calculated Qpp('^E) 
are zero (not shown), i.e. the electrons are on different 
dots, due to the Pauli principle. However, in large elec- 
trostatic quantum dot molecules^i^ (with Jee > Aee), 
where the lowest '^E state is expected to mix with higher 
energy triplet configurations and acquire non-zero double 
electron occupations, Q^^p > 0. (iv) QtotC'^u) is close to 
1 for all inter-dot distances (Fig. 0Jd) and both electrons 
are therefore occupying the same dots. 

As we have noted in the introduction, quantum com- 
puting based on two spins in dot molecules requires that 
(i) q['^1 ^1 i.e. the electrons should be localized on dif- 
ferent dots as much as possible. We now see this is satis- 



For 



fied for '^Eu, but not for the singlets ^E„ and ^Eg 

"'^Eg"'', we find Q[ot only at large inter-dot separa- 
tion d. (ii) Two spins should have high entanglement. 
Since entanglement is maximized for pure singlet and 
triplet states without double-occupancy (i.e. Qtot=0), 

this also requires that Q['^1 ^1. (iii) The inter-dot sepa- 
ration should be such that significant singlet-triplet split- 
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FIG. 4: The probability of two electrons occupying the same 
dot [Eq.lO] for (a) ^E^"', (b) ^E^, (c) ^E*''' states. Qtt 
(Qbb) is the probability of both electrons being on the top 
(bottom) dot, while Qtotai gives the total probability. 



ting exist. Conditions (i) and (ii) require large inter-dot 
separation d, while (iii) requires small d. Considering all 
the present results, the requirements for quantum com- 
putation are best met for the ground state -'^Eg"^ at an 
inter-dot distance of 6 to 8 nm, where we have significant 
Js-T, and Qtot is small. 

To conclude, we investigated the wave function asym- 
metry of both single-particle and many-particle wave 
functions in the quantum dot-molecules made of two 
identical, vertically stacked InAs/GaAs self- assembled 
quantum dots. For single-particle states, we find the 
asymmetry at short interdot separations introduced by 
the inhomogeneous strains. For many-particle states, we 
find that two electrons are always on different dots in the 
triplet state, while for the ground state singlet, the proba- 
bility of "bielectron localization" is not zero and increases 
rapidly as the two dots come close to each other. The 
ideal interdot separation where the two electron spins 
can be used for quantum information is around 6 -8 nm. 
At this distance the singlet-triplet splitting is large but 
the double electron occupation still low. 
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